GSI anomaly and spin-rotation coupling 
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We propose a model in which a recently reported modulation in the decay of the hydrogenlike 
ions 140 p r 58 + ] 142 Pm 60+ and 122 I 52 + arises from the coupling of rotation to the spin of electron and 
, nucleus. The model shows that the spin-spin coupling of electron and nucleus does not contribute 

to the modulation and predicts that the anomaly cannot be observed if the motion of the ions is 
rectilinear, or if the ions are stopped in a target. It also supports the notion that the modulation 
frequency is proportional to the inverse of the atomic mass and that no modulation is expected for 
the /? + -decay. 
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1. INTRODUCTION 

^ . Experiments carried out at the storage ring ESR of GSI in Darmstadt [1H3J reveal an oscillation in the orbital 
electron capture and subsequent decay of hydrogenlike 140 Pr 58+ , i42p m 60+ an( j 122^52+^ rji^g modulation has a 
period of 7.069(8) s, 7.10(22) s and 6.1s respectively in the laboratory frame and is superimposed on the expected 
' exponential decay. The "zero hypothesis" of a pure experimental decay has been excluded at the 99% C.L. and 
periodic instabilities in the storage ring and detection apparatus also seem improbable causes of the modulation. The 
effect has been extensively studied in literature 

We show, in the model proposed below, that a modulation arises in the probability that the system, initially in 
a superposition of hyperfine states (F = 3/2 and F = 1/2), finds itself again in such a superposition of hyperfine 
states after injection into the storage ring. The modulation has its origin in the spin-dependent part of the Thomas 
precession, and is compatible with the observed ESR modulation. The EC decay occurs for states with spin F = 1/2 
. because decay from the spin 3/2 state is forbidden by the conservation of the F quantum number [2|. We stress that 
the present paper differs in essential ways from Q because it takes into account all the relevant features of the GSI 
experiment, such as bound states kinematics, dragging effects, Thomas precessions of nucleus and electron and QED 
and derives the probability of the observed modulation from the time evolution of nucleus plus electron once this 
system is injected in the storage ring. 

The full Hamiltonian that describes the behavior of nucleus and bound electron in the external field B of the ring 
is H = Hq + Hi, where Hq contains all the usual standard terms (Coulomb potential, spin-orbit coupling, etc.), and 
Hi is (in units h = c = 1) 

Hi = -As ■ i - s • n e - 1 n n , (1) 
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where A ~ Z 3 ~ N x 10 14 Hz , N ~ 0(1) is the strength of spin-spin coupling, while 

3 

n e = u> ge +u>^ h -u>^ , (2) 

n„ = u gn + u%> - «w , (3) 

represent the precession of the electron spin and the usual spin precession of the nucleus in its motion in a storage ring 
that is assumed circular for simplicity. In and ([3]) , o; Sc n are the electron and nucleus spin precession frequencies 

due to the respective magnetic moments g e ^ n and Uc''"' are the angular cyclotron frequencies. The explicit expressions 
of all these quantities are given below. We refer ([l) to a frame rotating about the X3-axis in the clockwise direction 
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of the ions, with the 12-axis tangent to the ion orbit in the direction of its momentum and write B 
B = 1.197T is the GSI value (we are assuming that this is the average value over the circumference). 



Buz, where 



The Thomas precession iu^^' is related to the standard spin-rotation coupling that can be derived, for the electron, 
from the spin connection coefficients of the Dirac equation in a rotating frame 0, [To| . In our derivation, we neglect 
any stray electric fields and electric fields needed to stabilize the nucleus orbits, as well as all other effects which could 
affect the Thomas precession 11] . 

We indicate by and n the velocities of electron and nucleus relative to the lab frame. Using the composition 
of velocities, the Lorentz factor 7 = — 2 of the electron can be written in the form 7 = 7„7 e |„(l + II), where 

II = 0„ ■ 0„\„ = 0„0 r ,\„ cos 6*, 0„\„ is the velocity of the electron relative to the nucleus, i„ — 1/ yl — 0\ 
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0n ' 0e\n 
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0n0 e \n cos 6*, e i n is the velocity of the electron relative to the nucleus, 7„ 
. The explicit expression of is also useful 



n, and 



= 



1 



7n(i + n) 



e 



7gn 

7n + l 



0n + ln0n 



(4) 



The Thomas precession of the electron in the lab frame is given by u)Th 



7 2 d0 



A 0. The field B in the lab 



7 + 1 dt 

frame (where E = 0) is transformed to the nucleus rest frame and gives E' = j n n A B and B' = 7„B on account of 



/3„-B = 0. The equations of motion are 



d0e\n 



0e\n 0e\n ' ^e\n 



for the electron with respect to the nucleus and 



dPn 

dt 



dt n le\ n m 7e\n m 

/3„AB for the nucleus with respect to the lab frame. Here f e |„ = — e(E'+/3 e | n AB') = — e7„(/3 e |„+/3„)AB. 



Using dt = 7„(1 + U)dt n , taking e \ n • B = 0, E e |„ A e \ n = and E e |„ A j3 n — (averaged over the decay time of the 
ion in the storage ring), we find = — — — - — =r-[/3 e i„ + n ] A B. Neglecting spin-orbit coupling 1 , u)Th can 



be written as 



dt 
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where 



7»n 2 

7n + l 



- Y 



(i + n)[ 7e |„ 7n (i + n) + i] 



Iq 



7e|„7„(l+n) +1 



Y 



2 eW [ ln {2 - cos 2 6) - sin 2 9} v _ 0\ n 0l sin 2 , 
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3(7™ + 1) 



The coupling of the electron magnetic moment fi e = — -f — s with the magnetic field is described by 

2 m P 



H 9c — —fj, e ■ B" = n e 



B 



7 + 1 



0(0 -B) 



The Coulomb interaction also contributes to the Thomas precession. It generates the spin-orbit coupling term in the Hamiltonian of 
the electron Hq ~ |m~r (dV/dr)s ■ L. However s ■ L = + 1) — 1(1 + 1) — s(s + 1)] vanishes when j = s for 1 = and j = I rt 1/2 for 

£ ^ 0, therefore -ffc = in the ground state. Moreover, the effect of A in n = ^j(p — eA) is negligible in the present context. 
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Keeping only the quadratic term in (3 e \ n and using (3{(3 ■ B) = 



u 3 J 



7^(i + n) 2 



B, we obtain 



H 9c =Tfi e B, 



( 7 e %(/3 e | n ; u 3 ) 2 
7(7 + 1) 



(6) 



and from it ds/dt = i[H 9c , s] = u) ge A s = T/i e A B which yields 



2m e 



TB 



(7) 



In order to refer the spin precession to the particle orbit, the effective cyclotron frequency u) c must now be subtracted. 
Its value is obtained by computing the instantaneous acceleration df3/dt = uffi f}\L\. Omitting terms like a e |„ = 
^-E e |„, /3 e |„ A B and [B • (/3 e |„ A /3 n )]/3 n that vanish when averaged, as already pointed out, we find 



u> 



(e) _ 



3j B/3„l-(/3 e |„-ui) 2 QB f3 n „ 



m e P 7e|r 1 7n(l + n) 2 M 7 „/3 



u 3 , 



(8) 



where El 
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7n + l 



7»n 

7n + 1 



and 
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(1 + n) 2 



From ©,© and © we obtain 



7e|n7n 



where T is defined in (J6)) and 



1 - (gejn • AO 2 /3n „ 

" 7 e| n 7™(l + n) 2 /3 ' 



C7 



A, 



r 



M 7„ 



0(1 + n) 



7»n 

7n + l 



(9) 



(10) 



Notice that the standard result f2 e = — eBa e /m e , where a e = (\g e \ — 2)/2 is the electron magnetic moment anomaly, 
is recovered in the limit Q = 0. 

The calculation of g e -factors, based on bound state (BS) QED, can be carried out with accuracy even though, in 
our case, the expansion parameter is Za ~ 0.4. The BS-QED calculation gives fl2, 13| 



LtMEW + £g( a ) (aZ) 



(11) 



where C^(aZ) - 5 + ^(aZ) 2 + |(aZ) 4 . From (HQ) we obtain the values a e ~ -0.065122, a e ~ -0.0682112 and 
a e = 0.0505352 for 140 p r 58 + j i42p m 60+ an( j 122 j 52+ res p ec ti V ely. The addition of more expansion terms Q does not 
change these results appreciably. 



Consider now the nucleus with spin I. The terms of Q are u) 3n = g„/^jvB, where ^a? 



\± 

2m, 
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and give 
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M \ 2 Z 



(12) 
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PROBABILITY AND MODULATION 



Let \I,mi)i and \s,m s ) s be the eigenstates of the operators I and s. The total angular momentum operator is 
F = s + I. The angular momentum F assumes the values F = 3/2, 1/2, mp =3 / 2 = ±3/2, ±1/2, and tof=i/2 = ±1/2 
because 1=1, mj = ±1,0 and s = 1/2, m s = ±1/2. By making use of the raising and lowering operators 
F± = I± + s±, we construct the normalized and orthogonal states 

3 3 \ -11 n 1 h 



3 1 

2'2 )F 



11,0), 



I h 

2' 2' 



11,1), 
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2' 2 



2' 2 



>F = |1,-1>J 



I -I) 

2' 2 7 



1 S - 

1 h 



||l,0> r ||,|> s - 



2' 2 /s 



2, ,11, /I, , 1 1, 
3l 1 '- 1 >'2'2>' + Y3 |1 ' 0)j 2'-2>- 



The (6 x 6) matrix with elements {(j)i\H\\(j)j) has the eigenvalues 



where 



A 2l3 = - T — , A 5l6 - j T — + - r - 



A± = 9A 2 ± 4^fi e =f 4An n + 4(Sl e - fi„) 2 



and the corresponding eigenstates \i) (i = 1, . . . , 6) 

|1,4) = M , |2,5) 
^4± 



|3,6) 



1 + Al 



<h 



l + Bl 
1 



V / i+^±" 



05 , 



1 + A 2 + 



where 



and 



A± = 



B± 



9A - 2fl e + 2Q n ± 3^AT 
4V2(O e -fi„) 

04 + 2fi e - 2f2 n ± Sy 7 /^" 

4V2(fi e -fin) 
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In the limit A 3> ft e ,n we obtain 

\i)-\<Pi), (13) 

and 

Ai - A 2 = A 2 - A 3 = = " e + 2 "" , Ax-A 4 = -n e + 2fi n , (14) 

A 5 -A 6 = " e ~ 4 " n . (15) 
Notice that in these expressions the ^4-terms coming from the spin-spin coupling cancel out. 

2.1. Modulation induced by quantum beats 

These results must be now applied to the GSI experiment. Since the heavy nucleus decays via EC, only the states 
with F — 1/2 are relevant. For simplicity we confine ourselves to the Hilbert subspace spanned by the states {|5), |6)}. 
Here we follow [j| (see also @,0])- The decay processes involved in the GSI experiment 140 p r 58 + — > 140 Ce 58+ + v e , 
i42 Pm 60+ ^ i4d Nd 58+ + ^ and 122J52+ _^ 122^52+ + ^ can be scnematically represented as 

I ¥ + v e , (16) 

with obvious meaning of the symbols. At the initial instant t = (before injection into the ESR) the system nucleus- 
electron is produced in a superposition of the states {|5), |6)}, 



|I(0)> =X> a |a) = c s |5) + ce|6) 



with | C5 1 2 + |cg| 2 = 1. If one assumes, for simplicity, that the two states with energies A5 and Xq decay with the same 
rate T, at the time t the system evolves to the state 

|I(t)> = e- Tt ' 2 (c 5 e- jA5t |5> + c 6 e- lA6t |6)) . 

The probability of EC at time t reads 

PbcW =e- rt |(^ e ,F|5|I(i))| 2 = e- rt P EC [l + a 56 cos( W5 6t + 0] , (17) 

where (see Eq. (fl5|) ) 

^5 6 HA 5 -A 6 |= " e ~ 4 "" , (18) 

Pec = l( i/ e,I r |>S'|5)| 2 = K^FIS 1 ^)! 2 , a 56 = 2|cs| jes}, and finally S is the interaction operator 2 . The phase <r comes 
form possible phase differences of the amplitude ci and c 2 and of |(i/ e ,F|S|5)| and |(z/ e ,F|S|6)|. As (|17f and (|T8| show, 
the modulation of the decay probability does not depend on A. 



2 If we consider the Hilbert space spanned by the states 1131 . {|1), . . . , |6)}, then the probability 1171 1 assumes the form Pec{^) ~ 
[1 + 5Zi<j a ij cosaJijt], with LUij = \Xj — Aj |and i,j = 1,5. It contains 5 terms of which only one (that due to l|14[l and l|15|l ) contributes 
to the probability, while the others vanish because of EW selection rules. Assuming, therefore, that the states are equiprobable, the 
magnitudes take the value ay = -| ~ 0.2, with = 1(1,2), (1,3), (1,4), (2,3), (5,6)}. The values obtained in the GSI experiments 

are: o(Pr) = 0.18(3), a(Pm) = 0.23(4), a(I) = 0.22(2) [TJ|. 
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2.2. Estimate of 7 e 



We now compare the frequencies uj^/2tt, given by (fl8|) , with the experimental values ~ 0.14 Hz found for 140 p r 58 + 
and 142 p m 60 + an d ~ 0.16 Hz for 122 I 52 + and consider first the case uj^e = A5 — A6. We find 



\B 



3m,- 



T(a e 



le\nln 



u 



V z z 
.4. 



.4 



9n A 



A \ 2 Z 



2tt0.14 Hz. 



(19) 



where a bar on top means average values. These are computed by first expanding the quantities I e ,Q, T, U and V 
in terms of n < 1 and then averaging over the angle by means of (cos™ 6) = l 2( n +l) ■ Using (/3 e |n ■ u^) 2 = |/3 2 | n , 

i = 1,2,3, 7i Pr ' Pm ' I} = 1.43, 3i Pr,Pm) = 2.5, g$ = 0.94, up to 0OI 6 ) we obtain from (QU the numerical solutions 
(see Fig. [TJ 



7, 



(Pr) 

1 71 



1.07904, 7 



(Pm) 



1.08435, 7 



(i) 



1.05902. 



(20) 

(Pm) 



(Pr) 

which must be compared with the Lorentz factors of the bound electron in the Bohr model %\n ~ 1-0970, 7 e | n 

1.1040, and 7^ ~ 1.0776. The values (gDJ imply that the binding energies E = T + E p = -m[l - 7 e |„ + (aZ) 2 ], 
where T ad -E p are kinetic and potential energies of the bound electron, are given by 

£( pr ) 54.4keV, £ (Pm) - -58.2keV, E m ~ -46.3keV, 

in agreement with the values 



E^ r) = -49.5keV, E^ m) = -53.1keV, Jsjj? = -39.6keV, 



derived from the relativistic equation [2C 



En = 



RZ 2 



[aZf 



1-1- 
4n 



(21) 



(22) 



where R = 13.6057eV and n = 1. 



3. CONCLUSIONS 



In this paper we explain the GSI anomaly by means of a semiclassical model based on the Thomas precession of spins. 
The model has the following consequences: 1 ) It avoids all criticisms raised in fl6lj and in [Tt} because the Hamiltonians 
are essentially different. 2) There is no modulation in the /3 + -decay branch [18j. This because the Thomas precession, 
when computed only for a decaying charged nucleus gives rise to a frequency f2„ ~ eB/m p ~ 10 7 Hz and the 
probability ~ e~ r *[l + acos(O n t)]. The high frequency modulation term averages out to zero, and the probability 
obeys the standard exponential decay. 3) The GSI oscillations disappear when B = 0. 4) The model is consistent 
with experiments on EC decays of neutral atoms in solid environments that have shown no oscillations/modulations 
fl5j |. 5) The model predicts that uj^q ~ A^ 1 if the three terms on the l.h.s. of ([T9|) are of the same order of magnitude 
(notice however that the A-terms only affect the third digit of 7 e i„ and are not, therefore, a real discriminating feature 
of our model). 6) The model is consistent with the absence of a periodic transfer from active (F = 1/2) to sterile 
(F = 3/2) states. 
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